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The entropic uncertainty relation proven by Maassen and Uffink for arbitrary pairs of two 
observables is known to be non-optimal. Here, we call an uncertainty relation optimal, if the 
lower bound can be attained for any value of either of the corresponding uncertainties. In this 
work we establish optimal uncertainty relations by characterising the optimal lower bound in 
scenarios similar to the Maassen-Ufhnk type. We disprove a conjecture by Englert et al. and 
generalise various previous results. However, we are still far from a complete understanding 
and, based on numerical investigation and analytical results in small dimension, we present 
a number of conjectures. 


I. INTRODUCTION 

As a characteristic trait, quantum systems possess properties that are incompatible — properties 
that are equally real but mutually exclusive. In a pair of incompatible properties, if we have precise 
knowledge about one property, what we know about the other is necessarily imprecise. More 
generally, we can trade knowledge about one property for knowledge about the other and so know 
both imperfectly, and quantify the lack of knowledge by a suitable measure of uncertainty. Then, 
the compromises allowed by nature have their mathematical expressions in the form of uncertainty 
relations, which are inequalities that follow from the formalism of quantum theory. 

The study of uncertainty tradeoffs originated in Heisenberg’s pioneering work[T] of 1927 and 
was soon brought into a clear mathematical form by Kennard[2j. Weyl gave another early proof|3]. 
He was apparently unaware of Heisenberg’s paper and gives credit for the idea to Pauli, who seems 
to have learned it from Heisenberg in a letter prior to the publication of pQ . The modern textbook 
proof combining the Schwarz inequality with the commutation relations is due to Robertsonjl]. In 
this tradition the “uncertainty of observable X in the state p” is identified with the root of the 
variance of the probability distribution of the outcomes of an A-measurement on particles prepared 
according to p, i.e., 

,5A = y/tr(pA2)-tr(pA)2, (1) 

The key requirement for Heisenberg’s uncertainty relation 5Q6P > /i/2 to hold is that these 
variances are evaluated in the same state. The relation is thus a quantitative expression of the ob¬ 
servation that there are no dispersion-free states, and is hence of the type “preparation uncertainty 
relation”. This is in contrast to “measurement uncertainty relations” which express the feature of 
quantum mechanics that some observables may not be measured jointly, which also implies that 
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any measurement of one observable X implies a disturbance of the other in the sense that it cannot 
be inferred from a measurement on the state after an X-measurement. This aspect, although more 
prominent in Heisenberg’s paper than the preparation side, was made precise only recently[5] (also 

mil!). 

In this paper we are interested in preparation uncertainty relations for quantum systems of finite 
dimension d. A standard scenario in which this is of interest is the tradeoff between Welcher-Weg 
information and interference patterns at a multiport interferometer. In this minimalistic instance 
of wave-particle duality [8| one observable would detect particles on each of the internal paths of 
the interferometer, thus detecting a particle-like property, whereas the detectors at the end pick up 
wave-like interference. Uncertainty in this situation expresses the physical fact that if we prepare 
incoming particles so that they all go along the same path, we loose the interference contrast and, 
conversely, that large interference contrast is only possible when all paths are “traversed” with 
roughly equal probability. Another standard context for finite-dimensional uncertainty relations 
is quantum information theory, particularly quantum key distribution. Large parts of this theory 
have been developed in finite dimension, and there are many situations in which the incompati¬ 
bility as expressed by uncertainty relations plays an important role (e.g. in security proofs[9] of 
cryptographic protocols). 

What is common to these motivating instances of finite-dimensional uncertainty is that the 
outcomes of the respective observables are labelled in a completely arbitrary way. However, a 
variance depends not only on the abstract outcomes and their probabilities, but also on the real 
numbers we assign to them. For example, by multiplying all these numbers by the same factor we 
also multiply 6X. Moreover, variance will change if we permute the outcomes, which is as easy to 
do with beams in optical fibers as with freely re-codable bits of information. Basically motivated 
by such considerations, Deutsch[TU| suggested to use entropies to quantify the (lack of) sharpness 
of a probability distribution. This led to the famous entropic uncertainty relation established by 
Maassen and UffinkfTT]. to which we will refer to as the MU bound. It describes the sharpness 
tradeoff for the outcome distributions and py of two observables X, Y in the same state p in 
terms of their Renyi entropies Ha, Hp (see ([^), provided that these parameters satisfy the duality 
relation 
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When the observables X and Y are given in terms of their eigenbases {xt} and {yj}, so that 
p^(i) = {xi\p\xi) andp^(j) = {yj\p\yj), the MU bound is 

Ha{Px) + Hpipy) > -logmaxKyfclxj)!^ . (3) 


The bound becomes zero when the two bases share a vector, and maximal (namely logd) if the 
bases are mutually unbiased, so that all scalar products {yk\xj) have the same modulus. 

An alternative to entropies would again be variances, once one realizes that for defining a 
variance it is not really necessary to have M-valued random variables. It suffices to have outcomes 
in a metric space H with metric A, so that the variance of a probability measure /r on H becomes 

var(;u) = inf / p,{doj) A(a;,r/)^. (4) 

ri&n J 

When H = {1,..., d} the only permutation invariant metrics are A(i,j) = c(l — dij), and we will 
just fix the constant c = 1. Then 

var(p) = min > p{i) (1 — fc)^ = 1 — maxp(j). 
j j 

I 


( 5 ) 
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Up to a rescaling this is the so-called min entropy Hoo{p) = — logmaxjp(j). 

How then should we write an uncertainty relation in this general context? We will see that it is 
not wise to fix in advance the functional form of the tradeoff relation between Ha{px) and Hp[pY)- 
Instead, the best and most intuitive representation of the tradeoff is the diagram of all pairs 
{Ha{px)i Hi 3 {py)), ranging over all choices of input states p. An advantage of this representation 



FIG. 1: Entropy pairs for d = 2 and the observables X = and Y = [gx Y Left panel: The 

shaded set gives all pairs {H(p^), H(py))■ The contours describe the subsets which can be reached by pure 
states with a fixed admixture oi p = 1/2. Right panel: The shaded set is the “monotone closure” of the one 
on the left (see text). The solid curve represents the optimal bound: For entropy pairs on this bound it is 
impossible to reduce one entropy without enlarging the other. The thin line closer to the origin is the MU 
bound. 

is also that it changes in a simple way by a rescaling like the replacement of the variance ([^ 
by Hoc- For qubits (d = 2), all measures of sharpness are functions of each other, so all such 
diagrams are equivalent. Figureis drawn for the Shannon entropy H = Hi. Some details of the 
diagram of all pairs of entropies, shown on the left, are clearly not relevant for the uncertainty 
tradeoff, in which we ask how small we can simultaneously make the entropies. For this question 
it is the lower left corner of the diagram which matters, i.e., the set in the right diagram. It can 
be described as adding to any pair of entropies the full closed positive (north-east) quadrant above 
it. It is completely described by its lower left boundary, consisting of those entropy pairs with the 
property that for no other state one can have one entropy strictly smaller and the other at least 
as small. We consider the resulting curve as the complete description of the uncertainty tradeoffs 
between the entropies involved. Characterising this curve is the aim of this paper. 

We will always consider a quantum system in a d-dimensional Hilbert space, and consider two 
projection valued observables with d outcomes. This amounts to the choice of two bases {xi} 
and {uj}, and for the question at hand the choice is completely described by the unitary overlap 
matrix Uij = {xi\yj) modulo multiplication by diagonal unitary matrices or permutation matrices 
from either side. In the motivating standard case, closest to the case of position and momentum 
of continuous variables, the U represents the discrete Fourier transform of either the cyclic group 
of n elements or, if n is composite, another finite abelian group of order n. More generally, we 
also consider complex Hadamard matrices, i.e., unitary operators such that \Uij\ = Xj^fd for all 









4 


i,j. The bases are then called mutually unbiased, and we can think of a multiport interferometer 
generalizing a 50:50-beam splitter. Such bases also represent complementary pairs of measurements 
from the informational point of view. However, we will not restrict our study to these special classes 
of unitary matrices — several results will hold for arbitrary unitary matrices. For generalized 
observables (POVMs) or A:-tuples of observables similar questions can be asked, but we will not 
consider them in this paper. For the quantification of uncertainty or unsharpness we use the Renyi 
entropies Ha (1/2 < a < oo), and denote hy H = H^ the standard case of the Shannon entropy. 
Mostly we assume that the Renyi parameters a and (5 used for X and T, respectively, satisfy 
the duality relation Q. Again, the questions make sense also for other measures, e.g., related to 
majorization, or for variances, but these will not be considered here. We will also restrict ourselves 
to state-independent bounds, i.e., to the entropy pairs achievable by arbitrary states. When more is 
known about the state, for example about further expectation values, the entropy diagram for the 
subset may be quite different. Thus we do not consider inequalities like the Robertson inequality 
for variances, where the lower bound depends on the expectation of a commutator. 

Outline. In Sect. we briefly define all the relevant quantities and state our problem in precise 
mathematical terms. We present a brief review of previous results in Sect. |III[ In Sect. IV 


we 


provide a characterization of the case of equality in the MU bound and thereby show that the MU 
bound is not optimal in almost all cases. Our main results are presented in Sect. |V| We are not 
able to completely solve the problem in all its generality. However, we provide strong conjectures 
(Sect. VE) which, if true, heavily reduce the complexity of the problem. 


II. PRELIMINARIES AND NOTATION 


For a £ [ 5 , oo] the a-Renyi entropy of a probability distribution p G (0,1)“^ is defined by 


T^logEj=iP(jT if a 7 ^ 1,00 
Ha{p) = < - if a = 1 

— log maxj p{j) if a = 00 . 


( 6 ) 


The logarithms can be taken in any base (as long as it is always the same base). We follow 
the information theory convention of using base -2 logarithms, although base d would also be 
natural in our context, as it would normalize the range to 0 < Ha{p) < logd = 1. Monotone 
functions of the entropies tell the same story. In this sense we also cover “Tsallis entropies” 

r„(p) = (i-a)-i(i-E,p(jr)- 

Each entropy diagram will be drawn for a fixed choice of observables (i.e., bases) X,Y and 
values of the Renyi parameters a, /3, so that we consider a map / from the state space to given 
by 


f{p) = {h{p)j2{p)) = {Ha{p^x)^HpiPY)) ■ (7) 

Eor any choice we can define the order relation U on the state space, so that p Y p' stands for 
“/i(p) < flip') f 2 ip) < f 2 ip'y' ■ The uncertainty diagram is the monotone closure of the range 
{/(p)}, i.e., it is the set S containing precisely the pairs (/ii,/i 2 ) £ S for which there is a state p 
with fiip) < hi for i = 1,2 (compare EIG.[^. We call a state p optimal il p' Y p implies p C p', and 
hence f{p) = fip'). The corresponding optimal points in the entropy plane are characterized by 
the property that the uncertainty diagram contains no points to their south-west. We call the set 
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of all optimal points the curve of minimal entropies or the optimal bound. Therefore the optimal 
bound corresponds to a function 7 : ( 0 ,logd) —)■ ( 0 ,log(i) for which 

Haip'^x) (8) 

with the property that equality can be obtained for all possible values of Hp{py). 

If for some state the MU bound is saturated we call this state an equality state. The corre¬ 
sponding point in the entropy plane is an equality point. If an equality point exists we call the 
MU bound tight. The MU bound is said to be optimal, whenever it completely coincides with the 
optimal bound. 

A Hadamard matrix is a unitary matrix U with elements satisfying \Ujk\ = IjVd. The Fourier 
matrix is the matrix with components satisfying 

, j,k = 0,...,d-l . (9) 

The Fourier matrix is hence a special instance of a Hadamard matrix. This example generalizes 
to the wider setting of finite abelian groups, rather than just the cyclic group of d elements as in 
(§. To this end we consider the index set J for the first matrix index of U to equipped with a 
commutative binary operation turning it into a group. The second index is similarly labelled 
by the so-called dual group, denoted here by K. A symmetric way to express the relation between 
these groups is via the canonical bicharacter of the pair {J,K), which is a function : J x K ^ C. 
It has the property that the for every k the function j 1 —)• C(j)^) is a homomorphism from J to 
the complex numbers with modulus 1 , and that, conversely every such homomorphism is of this 
form for some unique k € K. Moreover, the same is true vice versa for the functions k 1 —>• Cifik) 
with fixed j G J. The Fourier matrix in this case is Ujk = d~^/‘^Q{j, k), where d = | J| = \K\. It is 
unitary and obviously a Hadamard matrix. When d is not a prime there are several non-isomorphic 
abelian groups of order d. 


III. PREVIOUS RESULTS 

There has been considerable work to generalize and improve the MU bound, e.g. by using more 
general entropy functions | 12 j or more than two observables HMS] (see also m for a review on 
entropic uncertainty relations). Most efforts, however, considered only the sum of the entropies 
(e.g. |18H25| 1. thereby already fixing the functional form of the tradeoff relation and not capturing 
all the information contained in the entropy diagram. 

In this work we are instead interested in characterising the curve of minimal entropies which 
we consider the optimal lower bound on the two entropies involved. There are, to the best of 
our knowledge, only very few results in the literature about the curve of minimal entropies in the 
finite-dimensional setting. In |26l |27j the authors note that the MU bound is not optimal in the 
simplest case of dimension d = 2 and compute the optimal bound for general unitary operators, 
but restricted to the Shannon case a = fi = 1. In[ 8 ]a conjecture about the entropy minimizing 
states is presented. We will see that this conjecture needs improvement. 


IV. EQUALITY IN THE MAASSEN-UFEINK UNCERTAINTY RELATION 

The MU bound provides a lower bound on the sum of two Renyi entropies that satisfy the 
duality relation Q. When characterising the curve of minimal entropies, it is natural to first 
investigate the case of equality in the MU bound. If the unitary operator linking the observables 
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is a Hadamard matrix, it is clear that the MU bound is tight. Indeed, any eigenvector of the 
observables, {xi} or {yi}, is an equality state. But can one also find equality points for arbitrary 
unitary operators? 

There already exist some results in the literature discussing this question, most importantly 
|28| and [29]. In the latter work the authors show the link between the two concepts of uncertainty 
principle and data processing inequality. Using this link the characterisation of all states that 
saturate the uncertainty relation reduces to the question of characterising all states for which the 
application of a certain channel does not imply loss of information. Employing this technique the 
authors can characterize all quantum states that saturate the MU bound in the restricted setting 
of observables related by Fourier transformation and Shannon entropies. A more general result 
was obtained in [28], namely a complete characterisation of all equality points in the special case 
a = /3 = 1, i.e. for Shannon entropies. Here we present an alternative proof of the uncertainty 
relation which allows us to generalize these from Shannon entropies to the case of arbitrary pairs 
of Renyi entropies that satisfy the duality relation. 

The main result of this section is the following Theorem. In its formulation the “support” of a 
probability distribution is the set of points with non-zero probability, and \M\ denotes the number 
of elements of a set M. 

Theorem IV.1. Let a,/3 > ^ be such that 1/a + 1//3 = 2, and let X,Y be bases with c = 
maxj^fc \ {yk\xj)\. Let p be a state, and denote by sx and sy the supports of the distributions pj^ 
andpy. Then equality in the MLf uncertainty relation 

Ha{Px) + Hp{p^) > log ^ (10) 

is reached if and only if p = |V')(V’I a pure state and, possibly after multiplying the basis vectors 
Xi,yj with suitable phases, the following condition holds: 

(xilV’) = \sx\~^^‘^ , (yj#) = and {yfxi) = c fori £ sx and j £ sy. (11) 

Moreover, 

, , , , 1 

sx Ry = ^. (12) 

Proof. We assume first that p = is pure, and will show that this choice is even necessary at 

the end of the proof. We fix '0 from now on, and choose phases for the basis elements so that, for 
i £ Sx, j & Sy we have 

= {xilif) > 0 and ipj = {yj\i>) > 0. (13) 

Note that this will change neither c nor the probability distributions. Furthermore, we assume 
without loss of generality that a < (3. We usually eliminate j3 by the duality relation, so the basic 
parameter to choose is a with 1/2 < a < 1. 

Our proof is inspired by interpolation theory, and involves the application of the maximum 
principle to a certain analytic “magic function” F. We do not pretend that finding this function 
is straightforward, since we also came by it in several stages of generalization and simplification. 
We define 

F{z) = c'-A' a,\y,) (14) 

i,j&s 


-1 


( 15 ) 
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where G s” is short hand for i € sx and 
j G sy, and (/?“ is the componentwise power of (/?, 
so that 


i 

and similarly for (p. The domain Q on which this 
function is analyzed is the strip 


^ = {z G C I 1 < Re z < 2}, 


(17) 


which is also depicted in FIG. Now since the 

sum (14) is finite and |r“^| = is bounded 

on Q for every r > 0, F is also bounded on Q, and 
the restriction of an entire analytic function. We 
claim that it is bounded in absolute value by 1. 
We estimate this separately for the two boundary 
lines. That is, for r G M we have, with Uij = {xi\yj) 




FIG. 2: Domain of F in the complex plane 


On the other hand. 


|F(l + ir)| = A 
= A 


ijSs 

= A ||^“||2 11^^112 = 1 . 


|F(2 + ir)| = c"^A2 




i,jGs 

hies 

< pf 

hies 


- 2/3 

i 


= A 


2\\,n\\2 


= 1 . 


(18) 


(19) 

( 20 ) 
( 21 ) 


Hence, by the maximum principle, |F(z)| < 1 for all z G Q. 

In order to relate this to entropies we consider the special value z = l/a, which always lies in 
the strip, but for a = 1 is a boundary point. We get 


= ipi 


ij&s 


j 

(22) 

J 

(23) 


(24) 
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where at (22) we used that <fiUij = (pj, and at (23) the definition of A and duality of a and f5. 
For taking the logarithm of this expression we use that 






and log(||^^||2^^) = 

and get, equivalently to F(l/a!) < 1, the inequality 

■ (log(c^) + < 0. 


, ^ , 1 \ 1 - a 

=-^ 


(25) 


(26) 


For a 7^ 1 we cancel the common factor and get the MU inequality. For a = 1 we always get 
F(l) = 1, and the MU inequality is obtained by taking the limit a —)> 1. However, it is better to 
express it instead by the derivative of F. For a = /3 = 1 we get 


F'(l) = -logc- < 0, 


(27) 


because for small e we must have F(1 + e) < 1. 

The advantage of this derivation of the MU inequality is that we have powerful characterizations 
of the equality case. So suppose that equality holds in the MU inequality. Then for a < 1 this 
means that F attains its maximum modulus 1 at the interior point 1/a of the strip Q, and the 
Phragmen-Lindelof Theorem |30j tells us that F = 1 is the constant function. For a = 1 we need a 
variant of the maximum principle due to Hopf[3T] (see, e.g. Thm. 2.7 in [32]), saying precisely that 
if the maximum is attained at the boundary with vanishing derivative we once again must have a 
constant function. In either case we conclude that F{z) = 1 for all z € G. 


With this information we can go back to the above estimates for (21), which must now be tight. 
The first step, the triangle inequality (19), is tight if all terms in the sum have the same argument, 
so up to a common phase the Uij for i £ sx and j £ sy must be positive. With the phase 
convention (13) this means Uij > 0 for all i,j in the supports. The second estimate ( |20[ ) is only 
tight when all Uij also have the maximum allowed modulus c. Hence Uij = c. If we consider U as 
an operator on vectors with support sy it thus maps to constant functions, so p must be constant 
on sx- By the same token (p must be constant on sy. Taking into account the normalizations we 
get all assertions of the theorem in the pure case p = |?/)('i/|. 

It remains to show that all equality states must be pure. So let V' now be any unit vector in 
the support of p and a = |V’)('!/|. Then we can write /? = Act + (1 — X)p' with A > 0, p' some 
other state, and similar convex relationships for the probability distributions. By concavity of the 
entropies, a must also be an equality state. Moreover, by strict concavity, a and p must have the 
same distributions and py = Pyi and hence the same supports sx, sy. Going through 

the proof for the pure equality state |V')(V'|, and in particular adopting the phase conventions made 
for ^ we find that Uij = c for all i £ sx and j £ sy. But then, if we apply U to any other unit 
vector ^|J' in the support of p we find that Uijj' is constant on its support sy. Hence ^|J' equals ^|J 
up to a phase, the support of p is one-dimensional, and p must be pure. 

An alternative proof of the necessity of purity, at least for the Shannon case a = /3 = 1, is via 
inequality [T2| 


1 


H{p''x) + H{p''y) > log ^ + H{p). 


(28) 


Clearly, for equality states the correction term, the von Neumann entropy H{p), has to vanish, i.e., 
the state must be pure. □ 
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An immediate consequence of Theorem IV. 1 is that for most overlap matrices no equality states 
exist, because 1/c^ is not an integer. Since the rows of a unitary matrix must be normalized, 
this integer is at most d, in which case we must have a Hadamard matrix. When 1/c^ < d one 
can build examples with equality by first solving a unitary matrix completion problem, starting 
from the known Sx x sy block. One then has to modify the matrix by unitary rotations on the 
complementary blocks so that all matrix elements become < c. The lowest-dimensional example 
is 2 = 1 !(? < d = 3, and the overlap matrix 


U = 


a a 0 
b —b a 
-b b a 


with a = and b = -. 

v/2 2 


(29) 


Some higher-dimensional examples can be generated by replacing the matrix elements a and b by 
aUi and bU 2 , where Ui, U 2 are any Hadamard matrices of the same dimension. 

By definition, Hadamard matrices have d orthogonal equality states with supports (|sx|, Isyl) = 
(1, d) and {d, 1), respectively. In prime dimension this is clearly the only possibility. However, even 
if the dimension is composite there may be no more than this, as the example j33| 



■ 1 

1 

1 

1 

1 

1 


1 

-1 

-T] 

—rf 

rf 

7 ] 

1 

1 


1 

rfi 

—rf 

rj2 


1 


rj-2 

-1 

jj2 



1 

jj-2 

—rj-^ 

j}-2 

1 

-rj 


_ 1 

jj-i 

rj-2 

_jy-2 

—rj-^ 

-1 


(30) 


with T] = ^ shows. Here one can mechanically check that none of the 300 3 x 2- 

submatrices has the property that all elements become equal after multiplication of rows and 
columns with suitable phases. Hence from Theorem IV. 1| it is clear that the point (log 3, log 2) on 
the MU-line is not accessible for any state. 

In the special case of a Fourier matrix (see the end of Sect.[T^for notations) we can get a complete 
description of the equality cases from Theorem IV. 1, as has been observed in Theorem 4.(1) of 


for the special case of a cyclic group. We will do the same for an arbitrary finite abelian group J. 
It turns out that the equality states are then directly linked to the subgroups of J and its dual K. 
The subgroups always come in pairs, i.e., when L C J is a subgroup, so is its annihilator|34j 


L^ = {k€K\yjeLC{j,k) = l}cK. 


(31) 


The basic result about annihilators is that (T-*-)-*- = L for every subgroup, so there is a ono-to-one 
correspondence between the subgroups of J and K, under which Li C T2 Lf- D L^- For any 
non-empty set L C J, we denote by xl the £^-normalized indicator function, i.e., xlU) = 
for j G L and xlU) = 0 otherwise. 


Corollary IV. 2. Let J be a finite abelian group, with Fourier matrix U, and L <Z J a subgroup. 
Then 


Uxl = Xl^, (32) 

and the vectors of the form x'(j') = C{j', k) xlU' — j), where j G JfL and k G K/L-^ are an 
orthonormal basis so that each |x^)(x^| is an equality state. Moreover, all equality states are of this 
form. 
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Note that in the formula for we can take arbitrary j € J and k & K, but two such choices 
(ji, ki) and (j 2 , k 2 ) define the same function when ji —j 2 G L and ki — k 2 G L-^- This observation 
is expressed by taking j, k in the respective quotients. 

We remark that, by the fundamental structure theorem of finite abelian groups, every such 
group is a cartesian product of cyclic groups, and has subgroups of every order which divides d 
(see Thm. 4.3 in [35]). Hence the equality points on the MU line are all points (log di, log ^ 2 ) with 
did2 = d. 

Proof. Let |V’)('0l be an equality state. The Theorem then says that for j G sx, and A: G sy we 
must have C{j,k) = fi{k)i'{j) for suitable phase-valued functions fj, : sy ^ C and u : sx ^ C. 
Now we can apply translations as in the construction of x' in the Corollary to get an equality state 
with 0 G sx and 0 G sy, from which we get ^{k)iy{fd) = 1 and /i(0)i^(j) = 1, so that the functions 
/i, V are actually constant. After applying an overall phase factor we can assume without loss of 
generality, that C,{j, k) = 1 for j € sx, and k € sy, and that if = Xsx- In terms of annihilators this 
is expressed equivalently by sy C or sx C Sy. 

When A: G we still have C{hk) = 1 for j G sx- But then {Uif){k) = {Uif){Q) > 0 and we 
must also have k £ sy. It follows that C sy. Combined with the already established reverse 
inclusion we get that sy = and, symmetrically sx = Sy- Note that since any set of the form A-^ 
is automatically a subgroup, we have shown that we can take sx = L, sy = L-*- for some subgroup 
L C J. 

We have so far only shown that Uxl is constant on L-*-, namely equal to y^|L|/| J|, coming from 
the summation of \L\ terms equal to and observing the overall normalization factor 



H(Py) 

log(6) 


log(3) 


log(2) 


log(2) log(3) 


log(6) 


H(Px) 


FIG. 3: Numerical sampling of the entropy diagram for dimensions d = 2 (light shading), d = 3 (medium 
shading) and d = 6 (dark shading) for Fourier-related observables and Shannon entropies. By Theorem 
|IV.1| the number of equality states corresponds to the number of divisors of the respective dimension. The 
optimal bounds (solid curves) are obtained by applying Conjecture V.6 and Conjecture V.7 presented in 
Sect.lVEl 
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of the Fourier matrix. We also have to show that C(j) k) = 0 whenever k ^ L-*-. However, in 

that case k induces a non-constant complex homomorphism on L, so it suffices to show that such 
functions add np to 0 on any finite abelian gronp. However, this is immediately obvious for cyclic 
groups, and hence follows for arbitrary groups by the structure theorem. So we conclude that Uxl 
is proportional to Xl-l, and since U is unitary, it must be equal, and |L-*-| |L| = | J|. 

Finally, let us count the translates x' for a given subgroup. Clearly, they are orthogonal to xl 
whenever either j-|-LnL = 0orA:-|- L-^ H L _L= 0. In other words, by taking one representative g 
from each class vuG/H and also one k from each class in K/L^ we get an orthogonal family. This 
has (I J|/|L|) (|iF|/|L-*-|) = | J|, i.e., is an orthonormal basis. □ 

For a product of abelian groups the Fourier matrix is the tensor prodnct of the Fourier matrices 
of the factors. Moreover one gets many equality states by tensoring, i.e., by taking subgronps of 
the form Li x L2 C Ji x J2. This additive structure is quite apparent from FIG.[^. It is therefore 
useful to note that this is also true without assuming the gronp structnre. This is shown by the 
following result. 

Corollary IV.3. Let Ui,U 2 be unitary operators of dimension di and d 2 , respectively. Suppose 
that for each unitary operator there exist an equality state and as characterized by Theorem 


Proof. First, note that max^^fc |(C/i (8> U 2 )jk\ = \Uijk\ maxj^fc \U 2 jk\- The MU relation then 

implies that, for any state cr on a di d2-dimensional Hilbert space, 


IV. 1 


Then the state cJeq = Ugq 1 


fjgq is an equality state for the unitary operator Ui0 U 2 . 


Ha{Px) + Hp{py) > -21ogmax|([/i ®U 2 )jk\ = -2 logmax max|t/2jA:| ■ (33) 

j,k j,k j,k 


In particular, for the state cJeq = dg, 


'eq 


(Tgq, we have 




2 

eq'i 

'X 


^eq\ 


^eq\ 


=-21ogmax|t/ijfc| max|t/2,jfc| 

j,k j,k 


(34) 


Hence, cieq is an equality state for Ui 0 U 2 . 


□ 


This Corollary should not be taken to suggest that only products will be equality states. For 
example, take the Fourier matrix of any abelian group of the form JxJ, which is the tensor product 
of two copies of the Fourier matrix of J. Then each subgroup L with | J\ elements generates a basis of 
equality states for the point (log | J|, log | J|). These are tensor product states for the subgroup L = 
{(j,0)|jG J} = J X {0}. But for H = {{j,j)\j G J} we get a maximally entangled eqnality state. 
Again, the basic idea of this example generalizes to more general settings. If Ui is any Hadamard 
matrix and Ui its complex conjugate, the maximally entangled vector Ylj \jj) is invariant 

nnder U = Ui 0 Ui. Hence both tp and U^p = ip belong to the equidistribution on d points, and 
\'ip){'tp\ is an equality state with entropies (logd, logd), just like \(p){<p\ with cp = Ylj |lj)- 
Perhaps one of the more surprising aspects of Theorem IV. 1| is that neither the characterization 
of the equality states nor indeed the value of the lower bound depends on a, pi. Hence we have 


Corollary IV.4. Let cJeq be an equality state, i.e. it saturates the uncertainty relation for some 
a, /3 > 2 satisfying the duality relation. Then aeq is also an equality state for all other pairs (a, pi) 
that satisfy the duality relation, including (a, pi) = (1/2,00), (00,1/2). 

The boundary cases for the ineqnality are proved by taking the limits on {a, pi), and since the 
lower bound is independent of these, equality carries over. However, additional states may then also 
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HpiPv) 



FIG. 4: Typical entropy diagram for Hadamard related observables in prime dimension for different values 
of a, (3 satisfying the duality relation ([^: a = 1/2 (light shading), a = 0.6 (medium shading) and a = 0.75 
(dark shading). The MU bound is optimal if and only if a = 1/2. 


satisfy equality. Indeed, Theorem IV. 1 does not hold in this case. As a counterexample consider 
an arbitrary Hadamard matrix U. Without loss of generality we can take it dephased, i.e., with 
all entries in the first row and column equal to Ij^/d. Consider then some arbitrary state r/ G 
with real and positive components to find 


max|(t/V>)fc|^ > = 2 


Taking the logarithm and using the definitions ([^ this is equivalent to 

logd > +Hoo{Py), 


(35) 


(36) 


which is > logd by the MU inequality. Hence all such states are equality states, and we can 

continuously interpolate between Hi = 0 and Hi = logd. Thus the MU bound coincides with the 

_ 2 2 _ 

optimal bound (see FIG.® and there is a continuum of equality states in contrast to Theorem IV. 1 


Another feature is true only in the boundary case, namely that for every U there is an equality 
state. To see this, let us consider an eigenstate Xj of X, for which Hi/2(7'x) = 0. But at the same 
time we have 


minHoo(py) = min(—logmax |(?/fc|xj)p) = —21ogc . (37) 

j j k 

One could summarize this by saying that in the boundary case {a,/?} = {1/2,oo} the MU bound 
is just too good to allow a useful characterization of equality. 
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/2 

FIG. 5: Consequences of concavity for the set of entropy pairs. 


V. CHARACTERISATION OF THE CURVE OF MINIMAL ENTROPY PAIRS 

Due to the study of equality in the previous section it is clear that the MU bound is, in almost all 
cases, not optimal, i.e. it does not coincide with the curve of minimal entropy pairs. To characterize 
this optimal bound is the aim of this section. We establish three general results that hold for 
arbitrary dimension: First, we prove that the curve of minimal entropies can be parametrized by 
pure states. Second, we show that for all real-valued unitary operators we can restrict the problem 
to real states. And last, we establish a necessary criterion for the Fourier case which all optimal 
states must satisfy thereby being able to characterize a whole class of potentially optimal states. 
Additionally, we provide a complete characterisation of the optimal bound for the simplest case 
of two-dimensional state space, d = 2. For d = 3 there is an analytic expression]^, which is well- 
confirmed by numerics, although not proved. However, for higher dimensions the optimal bound 
remains unknown. Nevertheless, we present random samples that suggest a number of conjectures, 
which, if true, vastly simplify the characterisation of the optimal bound. 


A. Sufficiency of pure states 

In this section we show that the optimal bound can be parametrized by pure states. At a first 
glance, this result may seem not too surprising since the situation is clear when minimizing only 
one concave functional f{p) over all states: In this case one can immediately restrict to pure states, 
since one of the convex components p' of p must always give a value f{p') < f{p)- However, the 
situation is not so simple when we consider a pair of concave functions, and the image of the state 
space under a two-component mapping / = (/i, / 2 ) as in Q. The direct consequence of concavity 
is then that for, say p = {pi + P 2 )/2, the point f{p) lies above the midpoint M = (^f{pi) + f{p2))/2 
in the coordinatewise ordering, i.e., fi{p) > (/j(pi) + /i(/ 02))/2 for i = 1,2 (see FIG. g). We 
therefore cannot conclude that the set {f{p)} is convex: the midpoint M is not in general in the 
set. Indeed this is clearly shown by the entropy diagrams, from which it is also clear that the 
complement is not convex either, except in simple cases. 

For the same reasons it is not obvious that it is sufficient to restrict to pure states. This is 
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highlighted by looking at the problem a bit more generally, considering the pairs of probability 
distributions in two bases. 

Proposition V.l. Consider two orthonormal bases X,Y in a Hilbert space and letp'^jPy denote 
the respective probability distributions in the state p. Then 

• If d = 2, then for every state p there is pure state a such that p^ = pf^ and py = Py- 

• If d < 3, then for every p we can find a convex decomposition p = pure states 

ai with p^ = p^fi for all i. 

For larger dimensions both statements fail. 

Thus, for d = 2 the range {/(/>)} is already exhausted by pure states, and for d = 3 the 
monotone closed uncertainty diagram can be computed just with pure states. For if f{p) is any 
point in the diagram, we can decompose into the a*, without any increase of /i, so by concavity we 
know one of the pure components has smaller / 2 . However, this proof strategy will fail for d > 4. 

Proof. (1) For d = 2, the set of quantum states p with the same distribution is the intersection 
of the Bloch ball with a hyperplane. Intersecting with the hyperplane for py we get a line, which 
also intersects the Bloch sphere, i.e., there is a pure state with the same distributions. 

(not 1) The example uses Fourier transform in d = 3. Two density operators have the same 
position distribution iff their diagonals coincide and the same momentum distribution iff the sums 
+ y) coincide for all y. Now consider a diagonal matrix with diagonal entries (1,1, 0)/2. 
A pure state with this diagonal will have just one non-zero phase in the 1-2 matrix element, so the 
sum with y = 1 will be non-zero other than for the mixed state. 

(2) Let us consider the convex subset K{p) of states with = p. We have to show that for 
d = 3 all extreme points of this set are, in fact, pure. Our method will also show that this fails for 
d > 4. 

First observe by just conjugating with a positive diagonal operator from right and left we get 
an isomorphism of K{p) and K{q), as long as p, q have the same support (of size d). So we may as 
well take p to be uniform, for which we write K{\) (Normalization factors are irrelevant here). 

Let us sort the potential extreme points by rank. Full rank is not possible, since then any vector 
with uniform distribution could be subtracted with a positive weight. Rank 1 is uninteresting, 
because it is of the form we want to exclude. This takes care of d = 2 and leaves only the rank 2 
case for d = 3. 

So let us consider the case of rank 2 for general d. Let <^i ,</>2 be two linearly independent 
vectors in the range of the density operator p = + \(l) 2 ) The condition that p has 

uniform position distribution means that |(/>i(x)p -|- |(/>i(x)p = 1 for all x. In other words, the 
pair <k(x) = (y(f)i{x), (l) 2 {x)) G is a unit vector for every x . Then we ask whether there is any 
non-zero vector T G of the form 'I'(x) = al(t>i{x) + a 2 (t> 2 {x) such that |T(x)| = 1 for all x. This 
would be a convex component of p with even distribution, so we could further decompose p. 

We can read this as a scalar product |(q;, <I>(a;))p. Think of the <l>(a:) and of a as represented 
on the Bloch sphere, where the geodesic distance is just a function of the above scalar product. So 
our question reduces to; Given d vectors on the sphere, can we find one further vector which has 
the same distance from each of them? 

Now for d = 2 this is obvious, and for d = 3 it works just like in the planar geometry of triangles: 
The locus of all points which have the same distance from $(1) and 4>(2) is a great circle bisecting 
their connecting geodesic at a right angle. Intersect with the bisector for ‘I>(2) and <k(3), which 
gives a point which has the same distance from all three points. Therefore, for d = 3, there are no 
extreme points of rank 2 , hence all are of rank 1 as claimed. 
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FIG. 6: States appearing in the proof of Theorem |V.2|as mapped to the two entropies plane. 


For higher d it is easy to find d points, which do not lie on a circle, i.e., there is no point 
equidistant from all of them. Hence there are extreme points of K{1) of rank 2. 

□ 

Surprisingly however, pure states can be shown to saturate all uncertainty diagrams, practically 
without assumptions on X,Y,a,j3. 

Theorem V.2. Let /i ,/2 he continuous concave functionals on the state space, define the order 
relation C as after equation 0 - Then for every state p there is a pure state a such that a T p. 

Proof. The plan of the proof is to show that for every non-pure p we can find another state a of 
strictly smaller rank such that a Q p. Then we can successively lower the rank, arriving finally at 
a pure state. 

Consider the face F of the state space generated by p. Its topological boundary dF consists 
precisely of the possible convex components of p of lower rank, and is connected. For each point 
a G dF there is a unique “antipode” . It is defined as 

0 -^ = (1 - A)cj) (38) 

for the smallest A for which the right hand side is positive semidefinite. It is clearly a state of 
reduced rank, i.e., G dF. We note that the required weight A cannot be 0 or I. 

We need not consider the case that cr E p, since otherwise we have found the desired element. 
Therefore, by exchanging the functions fi and /2 if necessary, we may assume that /i((t) > /i(p). 
We cannot also have fi{cr^) > fi{p)- Indeed, this would lead to the contradiction 

/i(p) > (I - A)/i((t) + A/i((T^) > flip). (39) 

Now consider a continuous curve [0,1] 9 t i—)• 7 (t) G dF connecting a and , i.e., such that 
7 ( 0 ) = a and 7 ( 1 ) = (see FIG.|^. Since fi was assumed to be continuous the previous argument 
shows that, for some t, /i( 7 (t)) = flip)- 

If /2(7(f)) < / 2 (p) we have found the desired element 7 ( 1 ) C P- The non-trivial case to 
consider is therefore f 2 {fiit)) > f 2 ip), or p C lit). Let A G (0,1) be the weight so that p = 








16 


(1 — + Xj{ty. Then by concavity, for i = 1,2, 


fi{p) > (1 - X)fi{-f{t)) + Xfi{-f{ty) 
> {I - X)fi{p) + Xfi{'y{ty) 
i-e., fiip) > fi{-f{ty). 


Therefore 7 (f)^ E p. 


(40) 

□ 


B. Sufficiency of real states for real unitary matrices 


From the previous section we know that for all unitary operators the complete optimal bound 
can be parametrized by pure states. Now we show that if the unitary matrix linking the two 
observables is real-valued, then we can further restrict the set of states for the complete optimal 
bound to the set of real-valued vectors. In this whole subsection we fix the Hilbert space to 
be C'^ with componentwise complex conjugation, so that the real vectors C C'^ are naturally 
embedded. 


Theorem V.3. Let /i ,/2 be continuous concave functionals on the state space and their inputs 
linked by a real unitary operator 17reai- -^Iso define the order relation C as after equation Then 
for every state p there is a pure and real state a such that a Q p. 


Proof. The idea of the proof is to employ again the proof technique of Theorem V.2, i.e. decompose 


a state in two states with the desired property (in this case, real states) and use the concavity 
property of the functions. 

Let V’ £ be a pure state. Since we are interested in a decomposition into real states, it is 
natural to consider the decomposition 


= y/Xv -|- iVl — Xw 


(41) 


where v,w are the normalized real and imaginary part of if and A = | Re('0)P ranges from 0 
to 1. We are only interested in the case where neither v ^ nor w ^ fi, otherwise the statement 
follows immediately. Furthermore, we assume without loss of generality that fi{v) > /i(V’)- Similar 


to the proof in Theorem V.2 we cannot also have that fi{w) > fi{4>) because we would then find 
the contradiction 


/i(V') > Xfi{v) + {l- X)fi{w) > fi{w) . (42) 

Consider now the states 

ip{t) := e'V (43) 

and their normalized real and imaginary part 

7 (t) :=Re{(p{t))/\Re{(p{t))\ , 

a{t) := Im(v5(t))/I Im((/j(t)) | (44) 


such that 


ip{t) = ^/pfij'yit) + iVl - p{t)a{t) 


(45) 







17 


where = || 7 (i)||. Note that fi{(f{t)) = fiiijj) for all t G (0, 27r). Also note that for a real-valued 
unitary operator the probability distributions and Py have the same form 

Px?Y = Pit)Px?Y + “ Pi^))Px?Y ■ ( 46 ) 

Due to continuity we know that there exists to such that either 7 (to) E V’) from which we obtain 
the desired statement, or ?/) C 7 (fo)- Using the concavity of the functions fi, the latter then implies 

fiW = > p(fo)/j(7(io)) + (1 - f^{to))fi{cr{to)) 

> /^(fo)/i(V’) + (1 - , (47) 

from which obtain fi[a{to)) < or equivalently cT(to) E V'- E 


C. Variatonal method 

So far we characterized the optimal bound by the order relation C. Equivalently, we may also 
consider an optimisation problem as mentioned in Q: Given some fixed value of 77^ (py) = 6 the 
optimal bound 7 is described by minimising Ha{Px)i i-O- 


7((5) = mm{i7„(p^)|i7/3(p^) = <5} 


(48) 


where 6 ranges from 0 to logd. However, performing this optimisation is in general quite difficult, 
especially because a nice characterisation of the constant entropy set {p|i7^(py) = <5} is not known. 
Instead, we restrict to optimising over a subset of this constant entropy set, namely states with 
varied phases. Clearly, this method will not yield a sufficient criterion for a state to be optimal. 
However, it provides us with a necessary criterion which allows us to identify a whole class of 
candidates of optimal states. 

More concretely, using Theorem V.2 we consider pure states p G and denote the components 


of the phase-varied state in Y basis by 


^3 = Pj exp ( -j- 0j 


27ri 


(49) 


for some phases dj. Varying these phases does not change the probability distribution, py = Py, 
and hence the phase varied states form a subset of the constant entropy set. For observables 
linked by Fourier transformation, we can optimize Ha{p\) over these states to find the following 
extremality criterion: 

Lemma V.4. Let the two observables X and Y be linked by the Fourier matrix (§ and let ijj 
denote an optimal state of this setup. Furthermore, let il) denote the Fourier transform off:. Then 
satisfies 


Im V'fc ^ 


dHa{p\)^ _ f2TTljk'- 


Proof. In order to optimize Ha{p%) we compute 

dHa{p 


fjj exp 


= 0 yk. 


7 ) 


dOk 


^ a77„(p^) 


6»=0 


i=i 


• |2 


d9k 


(50) 


= 0 


(51) 


6»=0 
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With w := exp (^) the Fourier transform of is defined as := Ylm=i hence, 


1 


d 


E 

m,n=l 




(52) 


Therefore we have 




dOk 


d 




0=0 m,n=l 

= ^ Im {ipkifjUJ^^) 


0=0 


Combining (51) and (53) we obtain the desired statement. 


(53) 

□ 


Any optimal state must necessarily satisfy the above criterion. This allows us to characterize a 
whole class of potentially optimal states: 

Lemma V.5. Let ip he a real-real symmetric state, i.e. a real state, p G satisfying the symmetry 
condition 


V>{j) = ^{d-j) Vj = l,...,d-1 


(54) 


or, equivalently, a real state with real Fourier transform, (p G 
criterion (50). 


Then p satisfies the extremality 


Proof. We first note a simple, but important property of real-real symmetric states: If is a real- 
real symmetric state and ^ is a state with components = f{pj), where / is any function taking 
real numbers to real numbers, then ^ is also a real-real symmetric state. For example, the Fourier 
transform of any real-real symmetric state is also real-real symmetric. 

Now p is assumed to be real-real symmetric. Hence, p is real-real symmetric. Define 




dHaipt) . 


(55) 


and note that ^ is also real-real symmetric. Importantly this implies that its Fourier transform, ^ 
is real. We therefore have for all k 


Im 


Pk 


dHa{Px 


d\pj 


2 Tjexp 


which finishes the proof. 


27rij k ' 


Im 



27rijA:^ 



(56) 

□ 


D. Simplest case: d = 2 

The results we presented so far are not sufficient to provide a complete characterisation of the 
curve of minimal entropy pairs. In what follows we therefore restrict to small dimension in order 
to reduce the complexity of the problem. 

More concretely, we investigate the simplest case, where the dimension of the Hilbert space is 
d = 2. In |26( [27] the authors characterized the curve of minimal entropy pairs for all unitary 
operators while restricting to the case of Shannon entropies. We now generalize their result to 
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FIG. 7: The optimal bound can be completely characterized in the qubit case (solid curves). The plot 
illustrates two entropy diagrams for randomly chosen unitary operators and entropy-pairs with a = /3 = 10 
(light shading) and a = /3 = 8 (dark shading). 


arbitrary pairs of Renyi entropies: First we show that for each 2x2 unitary oper ator U there 
is a real unitary operator U with the same entropy diagram. Then from Theorem V.3 we can 
immediately infer that the lower bound can be parametrized by real states. More concretely, 
our aim is to show that any unitary operator, which we can always write in {x,} basis up to an 
(irrelevant) global phase as 


/ cos{ip) sin((/ 9 )e \ 

Y — sin((/ 9 )e'® cos((/j) J ’ 


(57) 


is equivalent to the matrix 


/ cos((^) sin((/j) \ 
y — sin((/?) cos(<y 9 ) J 


(58) 


Indeed, the entropy diagram does not change if we first modify the unitary operator to U' = UV 
if R is a unitary operator satisfying Vxi = exp(i(/?i)xj for some phases tpi and all i, since then for 
any state p there exists a state p’ that yields the same pair of entropies. To see this, let p’ = pV 
to find that 

p^(i) = {xi\p'\xi) = {xi\V'^ pV\xi) = {xi\p\xi) =Pxii) (59) 

and 

Py'U) = {vjlp'ly'j) = {yj\vy^pyv^\yj) = {yj\p\yj) =pyU) ■ ( 60 ) 


u' = 


Now consider the modification 


cos{(p) sin((/?) 
sin((/?)e'^ cos((^)e^® 


(61) 
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obtained via the unitary operator 


^0 = 



(62) 


However, U’ yields exactly the same probability distribution as U. Hence, by Theorem V.3 the 
curve of minimal entropy pairs can be parametrized by real states. 

Since the real states form a one-parameter family it is not difficult to check that the states 


V’= (cos(0,sin(0) , (63) 

where the range of ^ is either (0, arccos(|[/i^i|) or (arccos(|C/i^i|), 7r/2) depending on whether 
arccos(|17i_i)| G (vr/d, 37r/4) or not, parametrize the curve of minimal entropy pairs for all uni¬ 
tary operators and all Renyi entropies. The problem is therefore completely solved in the simplest 
case d = 2 (see FIG. [^. 


E. Numerical sampling and conjectures 

In the previous section we characterized the optimal bound in the special case of dimension 
d = 2. To the best of our knowledge the problem is unsolved for all other dimensions. Instead 
the authors of [8] provide a conjecture stating that the curve of minimal entropies is traced out by 
states of the form 




(64) 


with pi-\- (d — l)p 2 = 1 in the case of complex Hadamard matrices and Shannon entropies. Due to 
the results of [26l[27] it is clear that this conjecture is correct for d = 2. The conjecture also holds 
true in the case d = 3 if we trust the numerics presented in FIG. where the solid curve directly 
corresponds to the states (64). However, for d = 4 we show that the conjecture already fails: For 


complex Hadamard matrices c = 1 /\/d and, hence, according to our analysis of equality in the MU 
bound there must be three distinct equality points, whereas the conjectured states only yield two 
equality points (see FIG.j^. 

However, we present two different conjectures which, if correct, explain how the bound in FIG.[^ 
and [H can be obtained: 


Conjecture V.6. (Product states for matrices with product form) 

Let the unitary operator U linking the two observables be a matrix of the form U = Ui®U 2 - Then 
for any state p there exists a product state pi ® p 2 with the same pair of entropies. 


The consequence of this our first conjecture is that the curve of minimal entropies for product 
form unitary operators in some composite dimension d = di d2 is just comprised of tensor prod¬ 
ucts of states that parametrize the curve in dimension di and d2, respectively. Indeed, from the 
additivity of the Renyi entropy it then directly follows that a state pd = Pd^ ® Pd 2 is optimal with 
respect to the unitary operator U = Ud^ 0 Ud 2 if and only if the marginals pd^ and are optimal 
with respect to the unitary operators Ud^^ and Ud 2 ■, respectively. We note that this conjecture also 


agrees with our findings for the equality states, especially Corollary IV.3 


Conjecture V.7. (Decomposition of the Fourier matrix) 

Let the two observables be linked by the Fourier matrix of composite dimension d = di d2. Then 
the entropy diagram does not change if we replace by <8) . 
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FIG. 8: Random sample of the entropy diagram for dimensions d = 2 (light shading), d = 4 (medium 
shading) and d = 8 (dark shading) for Fourier related observables and Shannon entropies. Our results 
falsify a previous conjecture by Englert et al. (dashed curves). Instead the optimal bounds are given by the 
solid curves, which are obtained by applying Conjecture |V.6| and Conjecture |V.7| 


The consequence of this second conjecture is that, although the Fourier matrix can, in general, 
not be decomposed into a tensor product of Fourier matrices of smaller dimension, the entropy 
diagram (and hence the curve of minimal entropy pairs) does not change under this replacement. 
Hence, if this conjecture were true, we could apply Conjecture V.6 and characterize the curve of 
minimal entropy pairs by states of product form, where the marginals parametrize the optimal 
bound in the respective smaller dimension. 

As an example let us consider Fourier related observables in dimension d = A. Employing both 
conjectures we know that it suffices to consider only the problem of characterising the optimal 
bound for Fourier related observables in dimension d = 2. But for such observables we already 
characterized the bound completely (see Sect. VD) and, hence, the optimal bound in d = 4 is 
traced out by product states with marginals given by (63). Indeed, this result agrees with the 
random sample (FIG.[^. In FIG. we also show other examples, where the numerics validate the 
two conjectures above. 

Note that the above conjectures are statements about the case of composite dimension, effec¬ 
tively stating that for a large class of unitary operators one only needs to solve the problem in 
prime dimension. The prime-dimensional case, however, still remains a hard problem. But we can 
provide two further conjectures that, if correct, vastly reduce the complexity of calculating the 
optimal bound in these instances: 


Conjecture V.8. (Independence of the optimal states of {a, ft)) 

If p is an optimal state for any unitary operator and any a, /? > ^ satisfying the duality relation 
([^, then p is also an optimal state for all other dual pairs. 

Note that we again excluded 


This conjecture can be seen as an extension of Gorollary |IV.4 
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FIG. 9: Random sample of the entropy diagram for real-real symmetric states in dimensions d = 2 (light 
shading), d = 4 (medium shading) and d = 8 (dark shading) for Fourier related observables and Shannon 
entropies. Restricting to real-real symmetric states does not yield the complete entropy diagram (grey), but 
seems to be sufficient to characterize the optimal bound. 


the extremal case {a,/3} = {1/2,00} for the same reasons as explained in Sect. IV, In FIG. the 
optimal bounds, although differently shaped, are traced out be the same states which supports 


Conjecture V.8 


The last conjecture only considers the case of observables linked by the Fourier matrix. 
Conjecture V.9. (Sufficiency of real-real symmetric states for Fourier) 

If p is an optimal state for the Fourier case, then there is a real-real symmetric state a as given by 


(54) with the same entropy pair. 


According to this conjecture it is sufficient to analyse the problem only for real-real symmetric 
states, which yields a huge simplification in both analytical and numerical treatments of the prob¬ 
lem. As an example consider Fourier related observables in dimension d = 3. If Conjecture V.9 
were correct, we already knew a characterisation of the optimal bound, since the real-real sym¬ 
metric states in this case form a one-parameter family and therefore trace out the desired curve. 
Indeed, for d = 3 the real-real symmetric states coincide with the states conjectured by [8] which, 
as mentioned above, trace out the bound if we trust numerics. FIG. also suggests the validity of 
Conjecture |V.9 

Furthermore, we note that real-real symmetric states are closed under the tensor product, in 
the sense that any tensor product of two real-real symmetric states is again a real-real symmetric 
state. Hence, Conjecture V.6|and Conjecture V.9 agree with each other. 
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VI. CONCLUSION AND OUTLOOK 

We investigated the curve of minimal entropies that completely describes the entropic uncer¬ 
tainty tradeoff between two observables. We showed that the lower bound on the sum of two 
entropies as given by the Maassen-Uffink uncertainty relation is not optimal in almost all cases 
and hence does not correspond to the curve of minimal entropies. To show this, we presented a 
novel proof of the MU bound that allowed us to analyse the case of equality in the uncertainty 
relation. 

In order to characterize the curve of minimal entropies, we provided three main results: First, 
we showed that the optimal bound can be traced out by pure states. Second, the optimal bound for 
real-valued unitary operators can be traced out by real-valued pure states. And last, we presented 
an extremality criterion, which any optimal state must satisfy. Numerical and analytical results 
for the case of small dimension suggest a number of conjectures that, if true, lead to a drastic 
reduction of the optimisation space. The optimal lower bound could then be computed. 
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